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Abstract 

We investigate the effect of topological disorder on a system of forced threshold elements, 
where each element is arranged on top of complex heterogeneous networks. Numerical results 
indicate that the response of the system to a weak signal can be amplified at an intermediate 
level of topological disorder, thus indicating the occurrence of topological-disorder-induced 
resonance. Using mean field method, we obtain an analytical understanding of the resonant 
phenomenon by deriving the effective potential of the system. Our findings might provide 
further insight into the role of network topology in signal amplification in biological networks. 
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Many social, biological, and physical phenomena have been well understood on the 
topology of complex networks. Wide studies of complex networks have shown that 
dynamics taking place in heterogeneous networks is dramatically different from that in 
homogenous networks, such as epidemic spreading, synchronization process, percolation. 
In view of this, it is necessary that a new type of topological disorder should be introduced. 
However, there is lack of unique definition of topological disorder at present. Alternatively, 
we propose a definition of topological disorder via measuring the disparity of degree 
distribution. We numerically and analytically demonstrate an interesting phenomenon, 
topological-disorder-induced resonance, by using coupling forced threshold elements on 
complex networks. Our results might have potential importance in understanding the role 
of network topology in the process of signal amplification of real biological networks, such 
as neural network. 
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Introduction 

It is well-known that many social, biological, and physical systems can be properly described by 
complex networks whose nodes represent dynamical individuals and links mimic the 
interactions among them [1]. The study of complex networks has indicated that many real-world 
networks exhibit small-world [2] and scale-free [3] features that are neither regular nor 
completely random. On the other hand, many scholars have laid their attentions to dynamical 
process happening on top of complex networks. Specially, vast researches in this field have 
shown that dynamical behavior on heterogeneous networks is dramatically different from that 
on homogeneous networks [4]. For instance, topological disorder could lead to a vanishing 
percolation threshold [5], the whole infection of disease with any small spreading rate [6], the 
Ising model to be ordered at all temperatures [7], the ttansition from order to disorder in voter 
models [8], synchronization to be suppressed in oscillator network [9], spatiotemporal chaos to 
be tamed [10], etc [11]. 

In such complex systems, there exists an important kind of disorder, namely topological 
disorder. Although a great many investigations have recognized the importance of this kind of 
disorder, there is lack of unique definition so far. In fact, one of the most major sources of 
topological disorder is the disparity of node degree, where node degree is the number of edges 
cormecting with the node. Therefore, it is feasible to characterize the topological disorder by 
standard deviation of node degree, as an alternative definition of the topological disorder. 
Disorder sometimes plays a counterintuitive role, for example, noise sometimes changes its role 
of conventional nuisance to a benefit. This is well-known as stochastic resonance (SR) [12], in 
which the right amount of noise is able to make a nonlinear dynamical system behave more 
regularly. Seminal works within the context of SR have been related to a weak forced bistable 
system [13] or excitable media [14] together with noise. When referring to the topological 
disorder, it is natural to ask whether the topological disorder induces a resonant behavior. To our 
best knowledge, it has not been reported previously in the literature. To the end, we here 
consider a system of coupled threshold elements, where the state of each element takes a binary 
value or 1, which is simultaneously decided by the states of its neighbors as well as an 
external signal. The system is modeled via networked structure where nodes denote threshold 
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elements and links connect them. We find that collective response of the system to the weak 
signal can be amplified at an intemiediate range of the topological disorder, thus indicating the 
occurrence of topological-disorder-induced resonance. 



The rest of this paper is organized as follows. In Section II, we propose the model and the 
approach to construct networks. Numerical simulation results of the model are presented in 
section III. Mean field analysis is put forward in section IV and main conclusions and discussion 
are addressed in section V. 



We consider a model of N coupled elements on underlying networks, where the states of nodes 
s. = 0, 1, ? e {1, • • • , N} are binary. The state is updated according to the following rule [15]: 



where w.p. denotes 'with probability', Q, is the set of neighbors of node i, and 0(D) is 

Heaviside step function. f(t) = A sin cot is the input signal, where A and CO are the amplitude 

and iiequency of the signal, respectively, h is the threshold and ^ is average degree. The model 
is a rather general paradigm for many real systems. The two states can be interpreted as being in 
favor an opinion or not, a neuron being firing or not, a gene being expressed or not, or several 
others. 

We first construct a network where node degree follows Gaussian distribution. The network is 
generated according to the MoUoy-Reed (MR) model [17]: each node is assigned a random 
number of stubs k that is drawn fi-om a specified degree distribution. Pairs of unlinked stubs are 
then randomly joined. This construction eliminates degree correlations between neighboring 
nodes. We ensure that the average degree K keeps constant and let the standard deviation of 
degree (J^ as the measure of topological disorder. Initially, each node is randomly assigned a 
state s(0)=0 or s(0)=l with equal probability 1/2. We perform numerical calculations by Monte 
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Carlo (MC) simulation. At each run, the first lO'^time steps are discarded to achieve steady state 
and the following lO"* time steps are used to investigate the system's dynamics. 

Numerical Simulation 

In Fig. 1 , we plot time evolution of the mean field m(t) = N ' .^^ .s^, (t) for different o-„ , with 

relevant parameters = 1 000 ,^ = 20,A = 0.5,^ = 0.28 and CO = 0.05 . One can notice 
that m oscillates with different fashions for different <7^ . For cr^ = or relatively small 

values of CTg , m oscillates around the value one (or zero when the initial value m{0) < 0.465 ) 

with the amplitude close to that of the input signal ^4. As cr^ increases, e.g. forcr^ = 6, m 

oscillates nearly in the whole allowable range [0, I] so as to dramatically increase the amplitude 
of the oscillations. With <J^ increasing again, m oscillates around a certain center value 

irrespective of m(0) , and both the amplitude and center value decreases as <J^ increases. 
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FIG. 1: Time evolution of mean field M for different <T^ on heterogeneous networks with 

Gaussian degree distribution. The other parameters are N = 1000 , K = 2.0 , h = 0.5 , 
^ = 0.28 and it* =0.05 

To quantify the response of the system to the input signal, we calculate the spectral 
amplification factor R , defined as the ratio of the output to input power at the corresponding 
driving frequency [16] 
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a =0, m(0)<0.465 
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m(t)e 
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where ((C^) and {\^ denote average over time and initial conditions, respectively. The 
dependence ofi?onO"^ for different y4 is shown in Fig. 2, with relevant parameters A'^ = 1000, 

K = 20,h = 0.5 ,A = 0.28 and O) — 0.05 . Here the amplitude of the input signal is set to be 
subthreshold, i.e. without topological disorder the response of the system to the signal is very 
faint. Each data is obtained via averaging over 50 different initial conditions and network 

realizations. With an increment of CJ^ , R reaches a maximum R^^ and then decreases, with R^. 

corresponding to a moderate magnitude of topological disorder, and thus indicating the 
occurrence of topological-disorder-induced resonance. As the magnitude of the signal A 
increases, the resonant peak becomes broader, the maximal height of the peak becomes lower, 
and the location of the peak slightly shifts to left. As the frequency of the signal 0) increases, 
R decreases monotonously and then vanishes when 6) becomes rather large, as shown in the 
inset of Fig. 2. 
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FIG. 2: i? as a flinction of (7^ on heterogeneous network with Gaussian degree distribution for 

A = 0.28 (square), A = 0.30 (circle) and A = 0.32 (triangle), with relevant parameters 
= 1 000 , K = 20 ,h = 0.5 , and O} = 0.05 . Symbol and solid line corresponds to simulation 
result and analytical one, respectively. The inset depicts i? as a function of CO with A = 0.28 



and CT^ — 5 fixed. 



5 



Mean Field Analysis 

In order to unveil the possible mechanism behind the above resonant phenomenon, we define 
as the probabilities that a node with degree A: is in state s = 1, and q as the probability that, 

for any node in the network, a randomly chosen nearest neighbor node is in state ^ = 1 . 
Furthermore, for any node the probability that a randomly chosen nearest neighbor node has 

degree k is A:/^ /K , where is degree distribution defined as the probability that a node 

chosen at random has degree k . It is supposed to be reasonable only in networks without 

degree correlation. The probabilities and q satisfy the relation 

q{t) = Y^kP,m,{t)lK. (3) 

k 

Note that q{t) differs, in general, from mif) = ^^^^^(0 ■ ^ particular, for all nodes are 

in state ^ = or in state ^ = 1 , one has q = m = and q = m = l, respectively. Suppose q 

is already known at a certain time t , one can calculate at the next time step. According to 
the evolution rule defined by Eq.l, in the absence of the external signal one has 

m,(t + V)=F(k,q)= X B(k,p,q), (4) 

p=[Kh] 

k\ 

where [□] is ceiling function, and B{k,p,q) = q'' (l — q) i% the binomial 

p\{k-p)\ 

distribution. Thus, we insert Eq.4 to the right-hand side of Eq.3, which yields the evolution 
equation of q 

q{t + \) = ^, {q{t)) = X kP,F(k, q(t))/K. (5) 

k 

The evolution equation of m is readily written as 

m(t + l) = {q(t)) = Y^P,F{k,q{t)). (6) 

k 

We now add the external signal to Eq.5 and Eq.6, which become 

q, m{t + 1) = (1- 1 fit) I) {qit)) + \ fit) \®{smcot). (7) 

We iterate Eq.7 and then calculate the spectral amplification factor /? as a function of (7^ , as 

indicated by line in Fig.2. It is clear that theoretical analysis well predicts the trends thati? 
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changes with and A . 

From time evolution of m (shown in Fig.l), there appears to be some clues which the system 
behaves bistable. Actually, by iterating Eq.6, we find that for initial value of the mean field 
/Mq > , m converges to one stable fixed point m^^,^ = , or else another one m^^ j ^ , 
which may imply the system is of bistable. Also, the evolution equation (5) has two stable fixed 
points and ^^^2 ' ^^^d an unstable one^^. In Fig.3(a), we show these fixed points as a 

function of cr^ . For any (7^ , one has m^^j = q^,^ = , mi^= q,^ and < ^,,2 • 

other hand, as (7^ is increased, , , 02^^,2 monotonously decrease, while ^,^2 decreases 

and then slightly increases. 
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FIG. 3: (a) ( q^ ), m^^,^ ( q^^.^ ) and m^,2 ( q„2 ) as a function of (T^ . (b) The effective potential 

of the system for three typical values of C7^ . (c) The height between potential wells and potential barrier 
AV (scheme in Fig.3(b)) as a function of CT^ . 

It is interesting to visualize the effective potential V{m) of the system, which can be done by 
converting the mean field equation to relaxational dynamics as 
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dfn dV (tn) 

= m{t-\-X) — m{t) = . For this purpose, we replace q{t) in Eq.6 with 

dt dm 

mif) + 1^(0 J where ^{f) = qif) — m{t) can be calculated by Eq.5 and Eq.6. Thus, the effective 
potential V{m) can be expressed as 

2 

V{m) ^- — r 2] (A:, iKh], m' + ^)dm'. (8) 

In Fig.3(b), the effective potentials for several typical values of are shown, with relevant 

parameters h = 0.5 and K = 20 . One can clearly distinguish that, for cr^ = 

(homogeneous network) there exist two asymmetric potential wells m^j , » separated by a 

potential barrier . If we apply an external signal to the system and guarantee the signal is 
weak enough, m will oscillates around one of the two potential wells (depending on the initial 
condition) with the amplitude comparable with that of the signal. As (7^ is increased from 

(Tg = , the potential barrier and the potential well locating at slightly shifts to left, as 
shown in Fig.3(a). Importantly, the potential barrier is lowered and the potential well locating at 
mw2 is raised, leading to the height between them, AFj and AFj , decrease (shown in 
Fig.3(c)). Thus the system may surpass it from one potential well to another such that the 
response of the system to the weak signal is amplified. When cr^ is further increased, A 

still decreases but AFj begins to increase. As a result, one potential well becomes very 

shallow but another one becomes very deep. With time evolution, the system quickly falls to the 
deeper potential well, and thus the system may fail to surpass the potential barrier again. 

It is useful to consider an altemative type of degree distribution, e.g. uniform degree distribution 
in which node degree is randomly selected in the range [A^ — A, + A] , where A is an integer 

between mdK — l. Similarly, the standard deviation cr^ =A/-\/3is the measurement of 
topological disorder. In Fig.4(a), we give simulation result and theoretical one of /? as a 



function of <J^^ for different A , which show excellent agreement between them. It is clearly 

shown that a maximum R arises at an intermediate value of <J^ . No matter which type of 

network we use, topological-disorder always gives rise to a resonant response. Moreover, we 
also find in this case the system has bistable potential. As shown in Fig. 4(b), we draw the 

bistable potential for three typical values of (7^, . Fig. 4(c) depicts the height between the 

potential barrier and the two potential wells, A Fj and A , as a function of cr^, . The 

obtained results show that the underlying phenomena and mechanism in 
uniform-degree-distributed network are similar to those in Gaussian-degree-distributed network. 
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FIG. 4: (a) 7? as a function of (T^, on heterogeneous network with uniform degree distribution for 
different A. Symbol and line correspond to simulation result and theoretical prediction, respectively, 
(b) The effective potential of the system for three typical values of CT^ . (c) The height between 

potential wells and potential baixier AI^ (scheme in Fig. 3(b)) as a function of (7^ . The other 
parameters are = 1000 , ^ ^ 20, A ^ 0.5 , and « ^ 0.05 . 



Consclusions and Discussion 

In summary, we show the influence of topological disorder on the collective response of 
coupled threshold elements to a weak signal on complex network, where topological disorder 
exhibits the disparity of node degree on the underlying network. We find that the collective 
response to the weak signal can be amplified at a moderate level of topological disorder. We 
analj^ically formulate mean field equation, and find the effective potential of the system is 
bistable in which the potential has two asymmetric potential wells separated by a potential 
barrier. As topological disorder is increased, the height of one potential well is always decreases, 
while the height of another one is first decreases and then increases. Only at a moderate level of 
topological disorder, the system surpasses the potential barrier such that the system makes 
amphfying response to the signal. Thus we gain good understanding of 
topological-disorder-induced resonant behavior. 

Topological disorder introduced in this paper stands for the disparity of node degree on the 
subject of complex networks, which can be easily distinguished itself from temporal disorder 
such as noise, and spatial disorder such as diversity in spatially extended systems [18, 19]. 
Although both noise [13, 14] and diversity [19] also induce a resonant behavior in bistable or 
excitable system, the starting point of our work is a different kind of disorder-topological 
disorder. Since the heterogeneity of degree is widespread in many biological network, our 
findings might give fiirther understanding toward the influence of heterogeneous topological 
structure in the ability to amplify weak external signals. 

Our numerical results and analytical treatment are based on the heterogeneous networks without 
degree-degree correlation. However, real-world networks often exhibit degree correlation 
feature [20]. Previous studies have shown that such degree mixing patterns have a significant 
influence in collective behavior of complex networks [20, 21]. Therefore, the effect of 
topological disorder in degree correlation networks on the present model deserves fiirther 
investigation. 
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